10 PRINT "CRITICAL DEPTH  COMPUTATION BY BISECTION METHOD": PRINT
20 READ Q, B, M, G: DATA 8.,2.,1.2,9.8:  PRINT "Q,B,M,G="; Q; B; M; G
30 DEF FNF (YY) = 1! - Q ^ 2 * (B + 2! * M * YY) / ((YY * (B + M * YY)) ^ 3! * G)
40 YMAX = 100!: YMIN = 0!: YCR = .5 * (YMAX + YMIN)
50 FOR J = 1 TO 14: X = FNF(YCR): IF X > 0! THEN YMAX = YCR ELSE YMIN = YCR
60 YCR = .5 * (YMAX + YMIN): NEXT J
70 PRINT "YCR="; YCR